Calculation of excited polaron states in the Holstein model 
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An exact diagonalization technique is used to investigate the low-lying excited polaron states in 
the Holstein model for the infinite one-dimensional lattice. For moderate values of the adiabatic 
ratio, a new and comprehensive picture, involving three excited (coherent) polaron bands below the 
phonon threshold, is obtained. The coherent contribution of the excited states to both the single- 
electron spectral density and the optical conductivity is evaluated and, due to the invariance of the 
Hamiltonian under the space inversion, the two are shown to contain complementary information 
about the single-electron system at zero temperature. The chosen method reveals the connection 
between the excited bands and the renormalized local phonon excitations of the adiabatic theory, 
as well as the regime of parameters for which the electron self-energy has notable non-local contri- 
butions. Finally, it is shown that the hybridization of two polaron states allows a simple description 
of the ground and first excited state in the crossover regime. 



I. INTRODUCTION 

Although more than half of a century has elapsed 
since Landau proposed that the charge carrier could 
be trapped by the distortion of a crystal lattice^ and 
Pekar introduced the term polaron^ a number of ques- 
tions regarding the single-polaron theory remain unan- 
swered. This is true even in the context of the Holstein 
Hamiltonian^ one of the simplest electron-lattice cou- 
pling models for the one-electron system. While the lit- 
erature pertaining to the ground state of such a model is 
extensive, much less attention has been paid to the ex- 
cited states, even at the low energies for which they are 
most interesting. 

The nature of the excited polaron states has been in- 
vestigated within the adiabatic approximation in Refs. 
I"i l,5l6l bv neglecting the polaron translation. These works 
provide a simple picture of the self-trapped polaron states 
for strong couplings. That is, the adiabatic softening 
of the phonon modes within the self-trapped polaron 
states results in several excitation energies below the bare 
phonon energy. It follows that the lowest excitations 
of the system are the polaron states for which the elec- 
tron and the phonons are strongly correlated. When the 
polaron translation is restored, each of the soft-phonon 
modes below the bare phonon energy can be expected 
to develop the corresponding band if the local dynamics 
remain adiabatic. Actually, moving polarons have been 
described within the adiabatic approximation by neglect- 
ing the force impeding the polaron translational motion 
due to the lattice discreteness^^ However, the band 
structure of the spectrum was not considered in these in- 
vestigations. For strong couplings, the bandwidth of the 
lowest band has been obtained by the adiabatic theory 
in the context of the simplest two-site modelii2*ii 

When considering the band structure it is important 
to realize that the polaron states are coherent in the 
range of energies below the phonon threshold (i.e., be- 
low the minimal energy for inelastic scattering) ^Sti^ One 
way to investigate their properties is by analyzing the 
single-electron and optical conductivity spectra. In the 



case of the single-polaron problem, the low-frequency co- 
herent part of the spectral weight directly determines 
the polaron energies. Most of the previous spectral 
investigations^^ 4 ^*^^^^ however, make no men- 
tion of the excited coherent polaron bands. The excep- 
tions in this respect are provided by those works em- 
ploying the dynamic mean-field theory (DMFT), which 
is exact in the infinite dimension limit. The DMFT re- 
sults predict one excited polaron band below the phonon 
thresholdi2i2i 

For the one-dimensional system, recent exact- 
diagonalization (ED) and variational approaches'^ ob- 
tained results for the lowest state of the first excited 
band. Provided that the local electron dynamics remain 
adiabatic, and that the adiabatic calculation of the self- 
trapped polaron energies gives several solutions below 
the inelastic phonon threshold, more than one excited 
polaron band is expected to occur in this energy range. 
The present paper shows that the ED approaches for the 
infinite lattice, as implemented in recent worksj22i"2£i 
can be extended to give very accurate results for the first 
few excited coherent polaron bands. In the range of pa- 
rameters for which the method converges, the coherent 
part of the spectrum is found to exhibit up to three ex- 
cited polaron bands, two more than previously described. 
These results provide a better understanding of the spa- 
tial and temporal structure of the polaron states, as well 
as of their symmetry (parity) properties. The latter can 
be used to rationalize the contributions of the bands to 
the low-energy single-electron and optical conductivity 
spectra, respectively. 

An additional and attractive feature of the excitation 
spectrum obtained herein is its potential to help clarify 
the physical picture of the polaron crossover for moder- 
ate values of the adiabatic ratio. It has been previously 
proposed^ in the context of the variational analysis 
that this crossover, in which the ground state evolves 
from a light to a heavy polaron state, can be under- 
stood as the anticrossing (hybridization) of two low- 
energy states. This question is reconsidered here in terms 
of practically exact eigenstates and their aforementioned 
parity properties. Specifically, the hybridization is found 
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to occur between states of equal parity, whereas states 
of opposite parity cross without any such hybridization. 
Furthermore, the hybridization of states in the crossover 
regime is characteristic not only of the Holstcin model un- 
der consideration. Recently, for example, on the basis of 
accurate numerical results obtained by the diagrammatic 
quantum Monte Carlo method^ for the Rashba-Pekar 
polaron model^ it has been suggested that the polaron 
crossover in that instance also involves a hybridization of 
several polaron solutions^ 

ED methods, analogous to the one employed here, have 
found widespread application in the ground-state calcu- 
lations of numerous other many-body problems. In this 
context, the present analysis of the excited states, to- 
gether with their symmetry attributes, is also of interest 
for the description of low-energy eigenstates and spectral 
properties from a more general point of view. 



increases. All states in the energy window are cor- 
related and will be referred to as polaron states. For 
E > £ (c) the electron and the additional phonons can 
form weakly-bound states, which results in a highly de- 
generate spectrum. 

The energies of the coherent polaron states are hence- 
forth denoted by E^, and the wave functions by l^jj- )■ 
K is the momentum of the polaron state (which is also 
the system momentum), while i denotes the band index. 
The lowest (ground) polaron band will be denoted by 
i = 0. In the present paper the term ground state is used 
for the state of minimal energy for a given momentum 
K, and not to the K — state only. Accordingly, the 
first excited state corresponds to the first excited state 
of a given momentum K, and so forth. 

For the polaron bands below the phonon threshold E^ 
the states • ) are the eigenstates of bq=o; 30 i 31 



II. HOLSTEIN POLARON PROBLEM 

The one-dimensional Holstein model of interest here is 
defined by the Hamiltonian 

H = -t^cj (Cn+l + Cn_i) + TlLJ^blbn 

n n 

- 9^ c « c « + M • (i) 

n 

t is nearest-neighbor hopping energy of the electron, huj 
is the energy of dispersionless optical phonons, while g 
is the electron-phonon coupling energy. c\ and 6^ are 
creation operators for the electron and phonon at lattice 
site n, respectively. As only the single-electron problem 
is treated, the spin indices have been omitted. 

The total momentum of the system (K) is the sum of 
electron and phonon momenta. As K commutes with the 
Hamiltonian , the complete set of states of the system 
can be constructed from eigenstates of H and K. The 
Hamiltonian is also invariant under space inversion. 
However, the analysis of the resulting parity properties of 
the eigenstates will be deferred until later in the present 
discussion [Sec. IIV D| . 

The eigenstates that will be considered herein are those 
falling in the energy window defined by 

£$L < E < £ (c) , £ (c) = E { °1 + Tiw . (2) 

E^ =0 is the minimal energy of the system (the zero- 
momentum polaron ground-state energy) and hw is the 
energy of the bare phonon excitation. When the electron- 
phonon coupling is absent (g = 0) the states below the 
phonon threshold E^ c ' are those of the free electron band. 
When g is switched on, the free electron states evolve 
into those states correlated with phonons. Furthermore, 
additional bands can appear below E^ as the coupling 
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More generally, it can be shown that for all eigenstates 
the simple sum- rule for the mean total lattice de- 
formation, given by 

X tot = $>„ = XoJ2«\( b n + b n)\K) = ^ , (3) 

n n 

is satisfied. In Eq. x is the space uncertainty of 
the free harmonic oscillator with frequency u>. Besides 
its physical meaning, the sum-rule for the mean total 
lattice deformation can also be used as a tool for checking 
the validity of results obtained with approximate polaron 
wave functions. 

III. EXACT TRANSLATIONAL METHOD 

In the case of the Holstein polaron problem, the ED 
approach uses only a finite number of states which con- 
tribute significantly to the exact polaron wave function 
for a given set of Hamiltonian parameters. By using the 
Hamiltonian matrix corresponding to the truncated (re- 
duced) basis, and the appropriate numerical scheme, one 
calculates the polaron wave functions and energies. The 
convergence of the results can be verified by increasing 
the number of basis states in the calculation. In most 
cases the results are very accurate, provided that the 
truncation procedure is well chosen. 

The ED method developed in the current paper is 
henceforth referred to as the exact translational method 
(eT method). The basis states of the eT method are 
given by«2£ 

|n ,n_i,m, ...,n m ) K = 
— ^2e iKja c^\n ,n- 1 ,n 1 ,...,n m ) j . (4) 
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The orthonormal wave function Q describes an electron 
which is surrounded by a cloud of phonons. The number 
of phonons at the mth lattice site away from the electron 
is given by n m , while K corresponds to the total system 
momentum and a is the lattice constant. The basis states 
i fi} of different momenta K are not mixed by the Holstcin 
Hamiltonian. 

If the adiabatic ratio t/Tito is not too large, the eT 
method can be used for studying polarons in the weak- 
and the strong-coupling regime, as well as in the crossover 
regime between them. The present paper is focused, par- 
ticularly, on the regime in which 1 < t/hcu < 5 and g is 
arbitrary. For g, t > Tito, the number of relevant basis 
states (@J becomes large and the problem of finding the 
polaron wave functions becomes untreatable. Unfortu- 
nately, the limit <?, t > Tiuj which is not considered here, 
is difficult for other known numerical methods as well. 

The matrix representation of the Holstein Hamiltonian 
in the eT basis leads to a sparse matrix. That is, it is 
straightforward to show that, by acting on the state (@J 
with the Hamiltonian Q , the maximum number of non- 
zero matrix elements per eT basis state is five. As an 
example, let us form the reduced Hilbert space of only 
five eT states, i.e., of the zero-phonon state |0)k, of three 
states with one phonon (phonon at the electron site, at 
the left site from the electron, and at the right site from 
the electron), 

|n = 1)k , |n_i = 1) K , K = , 

and finally, of the state with two phonons at the electron 
site | no = 2)^-. The corresponding Hamiltonian matrix 
is given by 
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In the limit g, t <C huj the eT method gives a good po- 
laron ground state with only these five states. Neverthe- 
less, for t ~ hui one usually has to work with a truncated 
basis of quite large dimension. In other words, the eT 
method generally requires a numerical scheme capable of 
dealing with large sparse matrices. 

For this purpose the well-known Lanczos algorithm 
appears to be the most appropriate choicei2Si2i Indeed, 
the previous papers using the eT method for ground- 
state calculations have employed this technique. An ad- 
ditional, attractive, feature of the Lanczos algorithm is 
that it is capable of finding not just one, but rather a 
few extreme eigenvalues (and eigenvectors) of sparse ma- 
trices, provided they are well separated (lying in the dis- 
crete part of the spectrum). Accordingly, the current re- 
sults have been obtained by the Lanczos procedure with 
the so-called local orthogonalization. In addition, the 



states have been calculated by the block-Lanczos proce- 
dure, with both variants giving the same results. 

It should be stressed that the present results are com- 
pared to the already cited ground-state results22i2i24 of 
the eT method. All states are checked through the sum- 
rule ©, as well as through their mutual orthogonality. 



IV. POLARON BANDS 
A. Numerical results 

At the beginning of our discussion it is instructive to 
observe how the polaron bands are formed in regard to 
the strength of the electron-phonon coupling g. For this 
purpose, the polaron bands are plotted in the four panels 
of Fig. as functions of the momentum (K ) , for t = 5 
and different values of g. As in the remainder of the 
paper, tuv — 1 is used as the energy unit. Notice that, 
as discussed in connection with Eq. J5J), only the part 
of the spectrum below the phonon threshold (between 

Ek=o an d E^) is shown. 

In the weak-coupling limit for the lowest band (-1 = 0), 
two regimes exist with respect to the critical momen- 
tum K C M£^ For K < K c the polaron state of energy 

E^' < E^ is the ground state of the system. For 
K > K c the ground state consists of the polaron and 
the unbound phonon excitation. The unbound phonon 
excitation carries the system momentum K, while the po- 
laron momentum is equal to 0. For K > K c the ground- 
state energy lies at the bottom of the incoherent part of 
the spectrum E > E^ c \ The first panel of Fig. ^ shows 
the eT weak-coupling results. Notice that the top flat 
part of the lowest band for K > K c has been cut by the 
frame. The reason is that the numerical error of the cor- 
responding eT states is slightly greater than that of the 
zero-momentum ground state, while the energy interval 
shown is exactly equal to the bare phonon energy fko. As 
the electron-phonon coupling increases, the lowest band 
becomes more renormalized, and finally, at some critical 
coupling, the whole band falls bellow E^ . 

By increasing electron-phonon coupling further, the 
additional polaron bands emerge below E^ . In the sec- 
ond and the third panel of Fig. the results for g = 2.8 
and g = 3.3 are shown, respectively. These are the 
choices of parameters that correspond to the crossover 
regime. For g = 2.8, only a part of the first excited band 
(i = 1) lies below E^ c \ whereas for g — 3.3 there are four 
polaron bands in the relevant energy window. From the 
third panel of Fig.^(g = 3.3) one sees that the top of the 
second excited band (i — 2) is at K — 0, and the bottom 
is at K — 7r I a. In addition, this band (i = 2) crosses the 
other excited bands (i = 1,3), i.e., the first excited band 
near the end of the Brillouin zone (for g < 3.3), and the 
third excited band near the center of the Brillouin zone. 

For strong couplings (g = 3.85) the eT results are 
shown in the last panel of Fig. ^ All the bands are 
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Figure 1: The polaron bands are plotted as functions of momentum (K) for t — 5 and four different values of g. Tiu = 1 defines 
the energy unit used throughout this paper. The results for g = 1 and g — 3.85 correspond to the weak- and the strong-coupling 
regimes, respectively, whereas results for g = 2.8 and g = 3.3 correspond to the crossover regime. Only the part of spectrum 
below the phonon threshold E < E^ is shown. 



very narrow. Although it is hard to distinguish the third 
excited band (i = 3) from the plot-frame, note that 

E m < E (c) for all K 

Figure [21 gives further insight into the polaron band 
formation. Here, the K = and K — n/a energy curves, 

shifted by E^L Q , are plotted for four bands (i is the band 
index). The results are presented for the crossover and 
the strong-coupling regimes. In the weak-coupling regime 
(not shown on Fig.[5J) becomes smaller than E^ 

for g w 1.6. 

It can be seen from Fig. [21 that the bandwidth of the 
lowest band decreases continuously with g. On the other 
hand, the bandwidth of the first excited state shows a 
more complicated behavior. Namely, for g « 3.15 this 
bandwidth is maximal, while, unlike for the lowest band, 
it decreases for smaller values of g. Its maximal value and 
the minimum of E^ =Q — E^ =a {g w 3.2) correspond to 
similar electron-phonon couplings. From Fig. [21 the band 
crossing which involves the second and the two other ex- 
cited bands can be clearly seen to occur in the crossover 
regime. That fact that these bands cross each other in- 
dicates that they belong to different symmetries, as will 



i=3 




Figure 2: The polaron bands below the phonon threshold 
shifted by E^ =0 are plotted as functions of g for t = 5 (hco — 



1). The band boundaries correspond to E^_ and E^_ n , a 
(except for the highest excited band for which this is only 
approximately true). 
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be discussed further below [Sec. IIV D| . 

In strong-coupling regime the spectrum consists of very 
narrow bands. Although in Fig. |2] the bands are well sep- 
arated, for very large couplings all of the excited bands 
approach from below. Finally, it can be seen from 
Fig. [3 that, for t = 5, the energy of the third excited 
band stays close to E^ both in the crossover and the 
strong-coupling regimes. The K values corresponding to 
the top of the third excited band depend sensitively on 
the Hamiltonian parameters. It should therefore be kept 
in mind that, in contrast to the other bands, the K = n/a 
curve only approximately determines the boundary of 
this band. In particular, the third excited band is slightly 
wider than suggested by Fig.|3]in the crossover regime. 



B. Strong coupling 

In order to understand the physical background of the 
bands shown in Figs. ^3 and [21 the strong-coupling and 
the crossover regime are examined separately in the fol- 
lowing sections. In the strong coupling regime the time 
scale relevant to the polaron translation is much slower 
than the time scale (1/w) involved in the local polaron 
dynamics. Therefore, the local interplay between the 
electron density and the lattice deformation is almost 
independent of the polaron momentum. From the en- 
ergy point of view the contributions corresponding to the 
polaron translation can be neglected, and the polarons 
can be treated, to a good approximation, as self-trapped. 
The deep potential well of the lattice deformation of the 
self-trapped polaron captures the electron. When the 
electron is light (hui < t), it is able to follow the slow 
motion of the lattice deformation (e.g., zero-point mo- 
tion), which results in the adiabatic renormalization of 
the phonon modes within the lattice deformation. 

The physical picture of the aforementioned self- 
trapped polaron follows from the adiabatic theory, and 
can be considered to be well understood. Thus, in 
order to achieve a better understanding of the cur- 
rent numerical results, it is convenient to compare the 
renormalized phonon energies obtained in the adiabatic 
limit to the spectrum calculated by the eT method. 
For the Holstein polaron problem the renormalized 
phonon modes have been calculated by different adia- 
batic approximations 4^*SiL2i2iii The procedure of Rcf. 
(Born-Oppenheimer approximation therein Sec. Ill) 
treats the lattice discreteness directly, while the polaron 
translation is neglected. This approximation is appro- 
priate for the present case as the lattice discreteness is 
important for the small self-trapped polarons, while the 
polaron translation has only a minor contribution to the 
energy. 

The three short-dashed curves in Fig. |31 are the renor- 
malized phonon energies of Ref. 0. The lowest excitation 
is a symmetric vibration of the lattice with respect to the 
polaron center (breathing mode). The next excitation is 
the antisymmetric vibration (pinning mode). The third 
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Figure 3: The energies of K — excited polaron states 
shifted by E^_ (solid curves) are compared to the energies of 
the renormalized phonon excitations obtained by the strong- 
coupling adiabatic approximation of Ref. (short-dashed 
curves). The long-dashed curve corresponds to E^L v i a and 
shows the lowest band narrowing. The inset displays only a 
small part of the spectrum below = 1. t = 2.5 (hoj — 1). 

excited mode is again a symmetric vibration, although 
extended over a larger number of lattice sites than in the 
case of the lowest excitation. In Fig. [31 the solid curves 
are the eT energies of the K — polaron excited states. 
In the strong coupling limit [the right part of Fig. [3] 
the results of Ref. are recovered asymptotically, from 
below. The inset of Fig. [31 where only a small part of 
the spectrum below E^ is shown, clearly demonstrates 
this behavior. For strong couplings, the positions of the 
excited bands are given by the the symmetric and an- 
tisymmetric phonon excitations of the adiabatic theory. 
In other words, the local dynamics of the self-trapped 
polarons is adiabatic. 

The energy Ejjr_ n j a ^ s plotted as the long-dashed curve 
in Fig. [3J This curve defines the bandwidth of the lowest 
polaron band which can be seen to become large in the 
crossover regime [the left part of Fig. [3], as the translation 
of the polaron becomes important. On the other hand, 
in the strong coupling limit, the eT method reproduces 
the narrow cosine polaron bands, 

E$-E%l n2t%[l-cos(Ka)}, 

where t^ ol is an effective polaron nearest-neighbor hop- 
ping energy for the ith. excited state. 

For strong couplings, the values of bandwidths may 
differ considerably from band to band. In Fig. Q] the 
bandwidths, as functions of g, are compared for constant 
A = g 2 /t fioj = 4.4. The solid curves are the results for 
the three lowest polaron bands. One sees that the ef- 
fective hopping t®, is increased if the polaron (i.e., the 
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Figure 4: The bandwidths (solid curves) of the three lowest 
polaron bands are plotted for 2 < g < 5 and A = g 2 /t Tilu = 
4.4. Each curve stops at the point where the top of the band 
crosses E^ c >. The bandwidths of the first and the second 
excited band are denoted by squares and circles, respectively. 
The long-dashed curve is the bandwidth of the nonadiabatic 
small polaron [a = 1 in Eq. JHJ]- In the inset, the cosine-like 
ground and two excited bands of the self-trapped polaron are 
shown, shifted by E^ =0 , as functions of momentum K (g — 4, 
t = 5, A = 3.2). All the energies are in units of hui. 

local lattice deformation) is excited. However, for strong 
couplings all bands are very narrow (4^ -C fuv), mean- 
ing that the polarons are very heavy (self-trapped). In 
real solids the coherent transport of such band states can 
be destroyed by the polaron-polaron interaction or by 
imperfections 

A close inspection of the eT results in Fig. 0] reveals 
that the hopping energy of the self-trapped polaron 
for the lowest band is given by the simple relation 

*S = * exph(.g/M 2 "(A)] . (5) 

Equation JSJ is an extension of the well-known expression 
for the nonadiabatic small-polaron hopping energy, which 
is obtained by setting a = liS. a(A) turns out to be a 
function of the adiabatic parameter A only, i.e., cv(A) is 
an adiabatic correction to the nonadiabatic small-polaron 
hopping. Notice that in the adiabatic limit A defines 
the spreading of the lattice deformation (the polaron size 
~ 1/ A) , as well as the renormalized phonon energies £ For 
A = 4.4 the excellent fit to the eT results in Fig. ^ (the 
short-dashed curve that follows almost exactly the eT 
results for the the lowest band) is achieved with a ss 0.86, 
even though the bandwidth of the lowest polaron band 
changes by several orders of magnitude. 

In summary, for the parameters under consideration 
[Fig. 21, it follows from Eq. (JSJ) that the translational dy- 
namics of the lowest band are essentially nonadiabatic, 
with adiabatic corrections. The narrow bandwidth de- 



fines the slow time scale for the polaron hopping. In 
contrast, Fig.|21shows that the fast local dynamics of the 
self-trapped polarons is adiabatic. The adiabatic soften- 
ing of the local phonon modes determines the positions 
of the ground and excited bands in the spectrum. The 
energy shift due to the softening is much larger than the 
bandwidths, which means that the polaron hopping can 
be neglected for the local dynamics. 

C. Crossover 

In spite of the fact that the almost exact results for the 
lowest polaron band in the crossover regime, obtained by 
various numerical methods, have been known for quite 
some time, the corresponding qualitative explanation of 
the polaron properties does not appear completely satis- 
factory. It has been suggested in Ref. |3{| by the use of 
global-local method^ that a simple empirical relation 
between Hamiltonian parameters 

g ST w huj + Vthu (6) 

describes the regime for which the fast change from light 
to heavy polaron ground state takes place. Indeed, gsr 
predicts, to a good approximation, the set of parameters 
for which the variation of the polaron effective mass (as 
a function of g) is the fastest. 

For t ^> hio the crossover occurs in the adiabatic regime 
when gsr ~ Vt Ucu, i.e., Xst ~ 1- Approaching Xst 
from the strong-coupling adiabatic side (A > Xst), the 
size of the polaron increases. Consequently, the Peierls- 
Nabarro (PN) barrier decreases and the tunneling of the 
adiabatic polaron to the neighboring sites becomes possi- 
ble. This effect is also responsible for the coupling of the 
pinning and the breathing lattice modest At A < Xst, 
the restoring force of the pinning mode due to the PN 
barrier can be neglected, and one may treat the polaron 
as a freely movingiiSiS Such a scenario thus describes 
the crossover from the self-trapped (heavy) to propagat- 
ing (light) polaron states in the adiabatic limit. On the 
other hand, for t <C huj, the nonadiabatic theory de- 
scribes the polaron translation for arbitrary g, and in 
particular, at g = gsT ~ ^ of Eq. ©. In this case, the 
polaron crossover appears essentially as a passage from 
the weak- to the strong-coupling ground state. 

A reinspection of Fig. [5] shows that the two well- 
separated time scales found in the strong-coupling 
regime, one related to the polaron translation (given by 
the bandwidth), and the other related to the polaron 
local dynamics (given by the renormalized phonon ener- 
gies), become comparable in the crossover regime. Fur- 
thermore, the behavior of the ground and excited bands 
for g « gsr indicates that the polaron translation plays 
an important role in the excitation spectrum. In other 
words, the influence of the polaron hopping to the neigh- 
boring sites on the renormalized modes cannot be ne- 
glected under these conditions, as it can for strong cou- 
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plings. The substantial difference in the crossover regime 
[see Fig. [3] between the eT results and the (adiabatically 
renormalized) phonon energies of Ref. is explained in 
this way. This difference, however, is not necessarily re- 
lated to nonadiabatic effects, as for the moderate val- 
ues of the adiabatic ratio under current consideration 
(1 < t/hui < 5) it is not clear to what extent the nonadi- 
abatic dynamics enter into the description of the polaron 
translation. 



D. Parity 

Irrespective of their temporal (or spatial) properties, 
when the polaron bands cross and/or anticross as in 
Fig. [21 their symmetry properties become important. As 
already mentioned, the Holstein Hamiltonian is invari- 
ant under the space inversion, and its eigenstates can 
be distinguished according their parity. A simple linear 
transformation relates the momentum and parity 

l^jf ) p eigenstates, 

\^ = \^)±P\^) = \^)±P\^%), 
so that 

P\*<$)=P\*%). (7) 

P denotes the space inversion operator, and P = ±1 are 
the even- and odd-parity eigenvalues, respectively. 

Considering the eT method, the parity can be di- 
rectly determined by inspection of the wave function 
properties. The eigenstates l^jf ) can be expanded in 
terms of the basis states (@J with expansion coefficients 

6kno, n_ i ,ni,...,n m , 

) = ^2 a «o,n-i,«i,...,n m |no, ...,n m ) K ■ 

Tli 

P acting on the basis state gives 

P |no, 71-1,711, :.,n m ) K = |n ,ni,n_i, ...,n_ m )^ K . 

Consequently, using Eq. Q, and the fact that 

is complex conjugate of | ) (time- reversal), one finds 
that the expansion coefficients satisfy 

0>n ,n-i,m,—,n m — P a no,ni,n_i ,...,n- m ' \°) 

a* denotes the complex conjugate of a. The expansion 
coefficients in Eq. (JHJ) stand for two local phonon config- 
urations, the first one is obtained from the second one 
when the phonons to the left and right of the electron 
are interchanged (n m —* u_ TO ). 



Applying the above analysis within the framework of 
the eT method yields P = — 1 for the second excited 
band (i = 2), while P = 1 for the ground and the other 
two lowest excited bands (i = 0,1,3). In this respect, 
the polaron bands inherit the symmetry of the renormal- 
ized phonon modes as obtained in the adiabatic limit. In 
Sec. lIVBl it was pointed out that the pinning mode cor- 
responding to the second excited band is an antisymmet- 
ric vibration of the lattice (having odd parity), whereas 
the other two lowest excited modes are symmetric vibra- 
tions (having even parity). Furthermore, one sees that 
the crossing (rather than the anticrossing) between the 
excited bands in the crossover regime, shown in Fig. |2 
involves bands of opposite symmetry under the space in- 
version. 

For K = and K = ir/a the linear transformation^ 
of the basis J2J 

|no,n_i,ni,...,ri m )jf = -y= (|n ,n_i,ni, ...,n m )K 

± |no,ni,n_i,...,n_ m )ii:) (9) 

defines two subspaces with different parities. It follows 
from Eq. JSJl that the P = 1 eigenstates I^k-o-k/o) 
long to the + subspace of Eq. (JHJ, while the P = — 1 
eigenstates belong to the — subspace. One sees that for 
K = and K = ir/a the parity actually defines the sym- 
metry of the local phonon configuration with respect to 
the electron. For P = 1 this configuration is symmetric, 
while it is antisymmetric for P = —I. 

E. Anticrossing 

In Ref. [2j| the correlated behavior of the ground and 
the first excited state in the crossover regime was ana- 
lyzed for moderate values of the adiabatic ratio using a 
method based on the variational approach. It was argued 
that the anticrossing of two physically different polaron 
states, one (heavy) for which the translation energy is 
almost negligible, and the other (light) for which this en- 
ergy is important, can describe in simple terms the mech- 
anism of the crossover. It follows from this interpretation 
that by increasing g in the critical region of parameters 
near g$T [Eq. JBJ], the contribution to the ground state 
from the light state decreases in favor of the heavy state 
which has lower energy at stronger couplings. On the 
other hand, the opposite change occurs for the first ex- 
cited state, which is heavier than the ground state for 
g < 9 st, while being lighter for g > gsT- 

In Fig. El the polaron ground and first excited energies 
obtained in Ref. |23j, denoted by CT, are compared to 
those of the eT method. Although the CT energies may 
be considered as a fair approximation, the anticrossing 
picture can be discussed more accurately in the context 
of the current eT results. One sees that the eT results 
confirm that the minimum 
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Figure 5: The ground- and first excited-state energies ob- 
tained by the CT method of Ref. I23I and the first excited- 
state energy obtained by the eT method, piotted as functions 
of g. The chosen set of parameters is the same as for Fig. 3 
of Ref.|H t = 5, gsr ~ 3.24 (hco = 1), K = 0. All curves are 
shifted by the ground-state eT energy. 



Figure 6: Difference of the energy, the mean number of 
phonons I1H . and the effective mass HUH of the ground and 
the first excited state for t = 1 (TkJ = 1). The minimum of 
E^L — E^P =0 and N^L = Nk =0 + 1 almost coincide with 
gsr = 2. Notice that for g < gsr the lowest excited state 
becomes heavier than the ground state (wig// > m i//)- 



states of the same symmetry (P = 1), corresponds to 
9 ~ 9 st- It is important to notice that the P = — 1 
states of the second excited band (i = 2) [Fig. [2] are not 
involved in the anticrossing. The rest of the investigation 
presented here is focussed on polaron effective mass and 
on the properties of the lattice deformation. 

(i) 

The polaron effective mass m 1 f can be evaluated nu- 



merically by using the relation 



eff 
18.22 



E 



(0 

K=Ak 



E 



(i) 

K=0 



(Ah a) 



(10) 



where Afc is a small finite deviation of the momentum 
from the value K — 0. m e i is the effective mass of the 
electron, and i denotes the band number. The ground- 
( TO i°f f ) an d the first excited-state (m^) f ) effective masses 



i (0) 

l eff 



> 



are compared for t = 1 in Fig. El Although 

n%ff f° r -9 ~ 9 st , m<, eff becomes smaller than nifh as g 
decreases. For larger t, the results are qualitatively the 
same [e.g., see Fig. [5], which means that the anticrossing 



picture indeed matches the behavior of m^fff ana - 
near gsT- 

For fixed g, the mean total lattice deformation (J3J is 
the same for the ground and for the excited states. How- 
ever, the local phonon cloud around the electron can be 
more or less localized, which affects the polaron hopping 
to the neighboring sites. A light state implies that the 
associated local lattice deformation is spread to a larger 
number of lattice sites than for the heavy state. Namely, 
such a deformation gives rise to a greater effective hop- 
ping integral. 



in 



(i) 
eff 



The mean number of phonons in the polaron state is 
given by 



K 



^blb m \* K 



(11) 



In the adiabatic limit, this number is quadratic in the 
local lattice deformation for the ground state. A more lo- 
calized (heavier) lattice deformation leads to larger . 
Within the adiabatic approximation the excited renor- 
malized phonon should, in general, increase the mean 

number of phonons with respect to . Particularly, in 
the strong-coupling limit one obtains 



hm N K = N K + 1 , 



(12) 



as the renormalization of the lattice vibrations becomes 
negligible. 

From Fig. H3 in which — is plotted as a func- 
tion of (7, one sees that the eT results tend to Eq. 1)12(1 for 
strong couplings. In the crossover regime, on the other 

hand, N% — deviates from Eq. i|12|) considerably. 

Nk is greater than N% + 1 for g < gsT and smaller 

than Nft + 1 for g > gsT- For larger i, the amplitude of 
the deviation increases even further. As shown in in the 
inset of Fig. □ N 1 k < for 3 - 25 < 9 < 3 - 75 > while 

the minimal and maximal values of N % — N ^ define an 
interval of almost 5 phonons. 
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-0.8 



0.6 



0.4 
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Figure 7: 



,. v(i,K) 

tion Xn 



The electron-lattice deformation correlation func- 
Eq. H13I . plotted as a function of g for n = 
and n = 1. In the former case X„ ,K ^ gives the correlation 
at the electron site, while in the latter it gives the correla- 
tion between the electron and the lattice deformation at the 
neighboring site. The difference between the mean number of 
phonons in the ground and the first excited state is shown in 
the inset, t = 5, K = 0, gsr ~ 3.24 (hu = 1). 



Additional insights into the local polaron properties 
can be provided by studying the appropriate electron- 
lattice correlation function, 



number of phonons, and the electron-lattice deformation 
correlation function) indicate that for t/tuo < 5 the an- 
ticrossing of two (light and heavy) polaron states occurs 



near gsr- Although criteria such as E^' =Q 
(o) r (o,o) _ Y (i,o) 



• ■ i (i) 
minimal, m eff 



EP =0 being 



o e ff, Xq ' — A ' , or alternative 
criteria, do not agree exactly, they all predict the polaron 
crossover to occur within the same very narrow parame- 
ter range, given approximately by Eq. ©. 



V. ELECTRON PROPERTIES 

A. Single-electron spectral function 

Denoting the vacuum state by |0), the zero- 
temperature single-electron Green function can be ex- 
pressed a a 41 ' 42 



G K (E) = (0\c K 



E - H + iQl 



Uo). 



For the energies below the phonon threshold E < E^ 
no polaron-phonon scattering takes place, which has the 
consequence that the imaginary part of the electron self- 
energy Y<k{E < E^) tends to zero Accordingly, the 
low-energy part of the spectral function Ak (E), denoted 
by Ak(E), is defined by the simple poles of Gk(E) at 



x (i,K) 



Cm (bl +n + b m+n )\^) . (13) 



X 



(i,K) 



is the mean lattice deformation induced at the nth 



site away from the electron. The correlation function 

is normalized in such a way that J2n X n K) = 1. As the 
correlation spreads to adjacent lattice sites, the correla- 
tion at the electron site Xq' K ^ decreases. Consequently, 

Xq' K ^ should be larger for the heavy than for the light 
polaron state. 

The ground- and first excited-state results for Xn' K ~ ^ 
are plotted in Fig.[7J The relationship between Ap 0,0 -* and 
Aq 1,0 ^ changes near gsr- This is exactly what one would 



for 



expect from the anticrossing picture. Aq 0,0 ^ < Aq 1 ' ^ 
9 ^5 9st- As the contribution of the heavy state to the 
ground state becomes dominant for g > gsr, the rela- 



tionship changes, and Ag 



(0,0) 



> A, 



(1,0) 



The spreading of 



the lattice deformation as a function of g may also be de- 
duced from x[ '' ^ and A^ 1 ' -*. The results corresponding 
to this lattice site, shown in Fig. [7] again lead one to the 
same conclusion. For g < g$T the lattice deformation 
seems to spread more for the ground than for the excited 
state, and vice- versa for g > gsr- 

All of the aforementioned ground- and first excited- 
state properties (the energy, the effective mass, the mean 




-lip 
E 

-12- 



-13 g=2.8 



-14 - 



-15 L 



M R =0 . 38 



g=3.4 
M^O.09 



g=3 



M =0 . 4 

g=3 . 2 



4 



A R (E) 



Figure 8: A^(E) given by Eq. fTH is plotted for t = 5, 
K = 0, and 4 values of g in the crossover regime. All spectra 
are broadened by a Lorentzian of width 0.05. For each curve 



the lower square denotes E = E 



(0) 



while the higher one 



denotes E = E (c) . Notice that for g > gsr ~ 3.24, the 
spectral weight associated with the first excited state is larger 
than the one associated with the ground state. Mk, given by 
Eq. 1151 , is the total normalized spectral weight of the ground 
and excited polaron states (below the phonon threshold). 
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E = E 



(i) 
K i 



A< (£ < = 



(14) 



i.e., the spectral function A^-(E) is denned by the po- 



laron energy E^S and the quasiparticle weight 



(i) 



K- \ C K/ 



\(9 

quasiparticle weight Z y ^J =Q , Q vanishes by parity for the 



In particular, for K — and K = n/a, the 



P 



-1 eigenstates |*^ = o, w / ) 

"if=0,7r/o' 



because the free electron 



(zero-phonon) state |ct r=n n/n ) belongs to the + subspace 



of Eq. ©. 

In Fig. \£\Ax =0 (E) is plotted as a function of E for a 
few different values of g in the crossover regime. The 
spectral weight corresponding to the ground polaron 
state is of the same order of magnitude as the one cor- 
responding to the excited (even-parity) polaron states. 
Moreover, Mr- defined by 



M K = J A<(E) dE , 



(15) 



1.5r 



0.2 



0.1 




3 g ST 3.5 



Figure 9: Z^_ and "/^ given by Eq. 1161 are plotted in the 
left and right panel, respectively (t = 5, gsr ~ 3.24, 7wj=1). 
The solid curves are the ground-state results, the long-dashed 
are the first excited-state results. In the left panel the short- 
dashed curves show the inverse effective mass m e i/m^Jj>. 



reveals that, in the crossover regime (for g < gsr), nearly 
40 percent of the total spectral weight confined to the 
ground and excited polaron states can be located in the 
energy window below the phonon threshold. Therefore, 
at least for g ~ g$T and for the values of adiabatic 
ratio under current considerations (t/huj < 5), the re- 
sults imply the existence of a few well-pronounced peaks 
in the single-electron spectral density below the phonon 
threshold. For stronger couplings, on the other hand, 
the quasiparticle weight below the phonon threshold is 
almost completely suppressed. 



B. Electron self-energy 

When the electron propagator Gk{E) is expressed in 
terms of the electron self-energy Hk(E), both Z^ and 
E$ of Eq. O can be related to £, K (E). E.g., the po- 

(i) 

laron effective mass m e ff is given in terms of the self- 
energy Yjk{E) by the standard formula 43 



/ W 
m e i/mlf f 



l + d eK E K (E) 
l-d E S K (E) 



z£Lo(l + 7 w ). (16) 



1 - 8e^k(E), and e K = K 2 aH 



where l/Z^ 

h 2 K 2 /2m e |. In Eq. l|16l) i is used to distinguish between 
results for the ground and excited (even-parity) states. 
7W is the abbreviation for the appropriately normalized 

(i) 

if-dependent contribution to m> e ff 

When the self-energy Y,x(E) is local (if-independcnt, 
0) the polaron effective mass is related solely to 
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(0 



(i) 

the quasiparticle weight Z K=0 . Thus, the non-local char- 
acter of the self-energy can be investigated by compar- 
ing 7W to unity^ Such analysis is particularly interest- 
ing in the context of DMFT, as the locality (implying 
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(') 



0) of the self-energy is an essential ingredient of 
this approach. 13 In both, the weak-coupling (g — > 0) and 
nonadiabatic (t — > 0) limits, T,jc{E) becomes localM^^ 
Within the Holstein-Lang-Firsov approximation! ^ and 
for t — ► 0, the self-energy Tijc(E) = £(-E) is, to a good 
approximation, given by its expansion to the first order 
in E over the whole energy range defined by the lowest 



band. In this case, Z 



(0) 
K 



Z^ defines the renormaliza- 



tion of the narrow cosine-like lowest band for any K. 
With the eT method can be obtained in terms of 



Ai) 



and 



l eff 



the latter being estimated numerically 



from the band dispersion at the center of the Brillouin 
zone [Eq. QlOfl]. The eT results agree with the afore- 
mentioned analytical findings for small t and/or gw^ It 
follows that significant non-local contributions to Y,k{E) 
can possibly occur in the regime where neither t nor g 
are small. In the left panel of Fig. |9] Z K=0 is plotted 
for moderate t as function of g, while the right panel 
shows 7W for the same set of parameters. Although 7W 
does not contribute to "m^ff substantially for g < gsr, at 
stronger couplings 7*- 1 -* reaches values of the order of unity 
for both the ground and the first excited state. That is, 
the non-local contribution to the electron self-energy is 
as equally important as the local one. However, it should 
be noticed that, in this regime, the quasiparticle weight 
zty becomes small. 
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C. Optical conductivity 

The interband optical conductivity may be evalu- 
ated within the linear-response theory from the current- 
current correlation function. The low-energy (coherent) 
contribution to the real part of the interband conductiv- 
ity at zero temperature is given byl^^l 

where are the K — polaron states, herein cal- 

culated by the eT method. The incoherent contribution 
to the interband conductivity denoted by a^(E), cor- 
responds to the continuum of the excitation spectrum 
above the phonon threshold. Since the current operator 
in Eq. JTJJ, defined by 

J = H ^ ] { c n+l c n ~ c JjC n +i) , 
n 

is odd under the space inversion, the nonvanishing matrix 
elements in Eq. H17|) are those between the P = 1 ground 
state |* (0) ) and P = -1 excited states |*W>. Actually, 
the incoherent part of the interband optical conductivity 
(Jft(E) involves only P = — 1 excitations as well. On 
the contrary, as argued in Eq. i|14|) , the spectral function 
Ak(E) contains information about the P = 1 part of 
the spectrum at if w 0. The optical conductivity at 
zero temperature is therefore described by those excited 
states for which the quasiparticle weight Zkzxo vanishes, 
i.e., those states that are not seen in the single-electron 
spectral function. 

The real part of the optical conductivity cjr(E) in- 
cludes the interband part and the (intraband) Drudc 
term at E « 0. It follows from the well-known sum-rule^ 
that the total spectral weight of <Jr{E) is given by the 
mean value of the electron kinetic energy in the ground 
state. Consequently, not only is the eT method capa- 
ble of calculating (E) , but the total spectral weight of 
<7r(E) is also accessible. Furthermore, one may estimate 
the weight of the Drude term from the polaron effective 
massi** 

According to the eT results, for t < 5 there is only one 
interband transition \'^ < ^ > =Q ) + — > I^k'Lq)^ which con- 
tributes at zero temperature to cr^(E) [Eq. l(T7l) ]. It is 
interesting to find the corresponding non-Drude spectral 
weight below the phonon threshold. However, E K ' =0 is 
less than E^ c > for relatively strong couplings, for which 
most of the spectral weight of ctr(E) belongs to high en- 
ergies. In this parameter regime only a few percent (or 
less) of the total spectral weight of (Jr{E) corresponds to 
<7^ (E), or to the Drude term. This simple low-energy pic- 
ture of ctr{E) becomes certainly more complex at finite 
temperature, particularly in the crossover regime where 
the bandwidths of the coherent polaron states increase. 



Studies at finite temperature, however, require different 
methods than the one presented here. 



VI. SUMMARY 

This paper reports an exact-diagonalization study 
of the ground and excited polaron states in the one- 
dimensional Holstein model. For values of the adiabatic 
ratio t/Tiu) < 5, accurate energies and wave functions are 
obtained for translationally invariant solutions of the in- 
finite lattice problem. The chosen method is restricted 
to the part of the spectrum below the phonon threshold, 
for which there are no phonon excitations uncorrelated 
to the polaron. 

The eigenstates of the Holstein Hamiltonian can be dis- 
tinguished according to their parity. This property fol- 
lows form the Hamiltonian which is invariant under the 
space inversion. It is shown that only the odd excited 
states of the system are relevant to the optical conduc- 
tivity at zero temperature. However, the contribution of 
the coherent excited polaron state is rather small with 
respect to the total spectral weight. On the other hand, 
it is the even states that contribute to the single-electron 
spectral function for if a: 0. In this case, the contribu- 
tion of the coherent excited polaron states is found to be 
important, particularly in the crossover regime. 

In the strong-coupling regime the results agree with 
the picture of self-trapped polarons, which may be re- 
garded as well understood. The spectrum is character- 
ized by very small polaron bandwidths. Consequently, 
the time scale associated with the polaron translation is 
very large, and the contributions to the fast local po- 
laron dynamics from the polaron hopping are negligi- 
ble. By making a comparison with the eT results, it is 
shown that the adiabatic theory (the Born-Oppcnheimer 
approximation) provides a good description of the local 
properties of the self-trapped polaron. Namely, the ener- 
gies of the renormalizcd phonon excitations are close to 
the eT excited-band energies. Two of the excited bands 
correspond to the symmetric adiabatic vibrations of the 
lattice, while one corresponds to the antisymmetric (pin- 
ning) vibration. Finally, it is found that for strong cou- 
plings the electron self-energy shows significant non-local 
(if-dependent) behavior at moderate t. 

As the electron-phonon coupling decreases at moder- 
ate t, the simple band structure found in the strong- 
coupling regime (i.e., the narrow and well separated po- 
laron bands) evolves notably. The effective hopping in- 
tegral of the ground and excited polaron states defines 
a time scale which, in the crossover regime, is compa- 
rable to the time scale relevant for the local polaron 
dynamics. In the crossover regime the excited (odd- 
parity) band which corresponds to the pinning vibra- 
tion in the strong-coupling regime, crosses the other two 
(even-parity) bands which correspond to symmetric vi- 
brations. The results suggest that the same critical set 
of parameters defined by gsr [Eq. ©], found for the po- 



12 



laron ground-state crossover, may also be associated with 
the rapid change of the low-energy excitation spectrum. 

Aside from the fact that the polaron translation and 
the local dynamics mix in the crossover regime, the in- 
termediate (moderate) values of the adiabatic ratio em- 
ployed herein present additional difficulties for qualita- 
tive understanding. Namely, for 1 < t/ftu) < 5, both the 
adiabatic and the nonadiabatic contributions could be 
important for the polaron crossover. Nevertheless, in the 
context of the eT approach, the analyses of the energy, 
the effective mass, and the lattice deformation properties, 
all agree in one important respect. That is, for moderate 
values of the adiabatic ratio, two polaron states of even 
parity (one heavy and one light) undergo an anticrossing 



in the same region of parameters for which the ground 
state crosses over from the weak- to the strong-coupling 
regime. 
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